中 国生 态 农业 学 报 2017 年 5 月 第 25 肉 第 5 期 RE 
Chinese Journal of Eco-Agriculture, May 2017, 25(5): 656-664 a 


DOI: 10.13930/j.cnki.cjea.160843 


任 志 远 ， 孙 艺 杰 ， 吴 林 筱 . 1990 一 2012 年 榆林 市 土地 生态 风险 综合 评价 四. 中 国生 态 农业 学 报 , 2017, 25(5): 656-664 
Ren ZY Sun Y J, Wu L X. Comprehensive evaluation of ecological risk of land during 1990-2012 in Yulin City[J]. Chinese 
Journal of Eco-Agriculture, 2017, 25(5): 656—664 


1990 一 2012 年 榆林 市 土地 生态 风险 综合 评价 ” 


任 志 远 , 孙 艺 术 , 吴 林 徐 
(陕西 师范 大 学 旅游 与 环境 学 院 /陕西 师范 大 学 地 理学 国家 级 实验 教学 示范 中 心 ”西安 710119) 


摘 要 : 土地 利用 / 履 盖 变化 对 区 域 生 态 安 全 和 维持 生态 系统 服务 功能 起 着 决定 性 作用 ， 土 地 生态 风险 评价 是 
衡量 土地 资源 可 持续 利用 的 重要 指标 。 本 文 以 生态 脆弱 的 陕西 省 榆林 市 为 研究 对 象 ， 依 据 榆林 市 土地 系统 的 
基本 特征 ,选取 表征 干旱 风险 的 降水 距 平 百分率 指数 、 表 征 人 类 活动 类 型 及 强度 的 土地 利用 结构 风险 及 表征 
土壤 侵蚀 的 土壤 水 力 侵 蚀 、 土 壤 风 力 侵蚀 4 个 风险 指标 , 构建 土地 生态 风险 综合 评价 模型 ,对 影响 榆林 市 
1990 一 2010 年 土地 生态 风险 的 因素 进行 定性 和 定量 的 评价 与 分 析 。 研究 结果 表明 : (1)1990 一 2012 年 榆林 市 的 
降水 距 平 百分率 指数 均值 哇 先 减 后 增 的 趋势 ,其 低 值 区 主要 分 布 在 偏 东 部 地 区 ; (2)21 年 间 土 地 利用 结构 生态 
风险 整体 上 呈 先 增 后 减 的 趋势 ， 区域 差异 较 大 且 变 化 明显 ; (3) 土 壤 水 力 侵 蚀 较 剧烈 的 区 域 主要 分 布 在 东部 及 
南部 地 区 , 21 年 间 均 值 呈 先 增 后 减 的 趋势 ; (4) 土 壤 风 力 侵蚀 较 严 重 的 区 域 主 要 分 布 在 北部 地 区 ， 随 着 时 间 变 化 
其 均值 一 直 呈 阅 减 趋势 ; (5) 土 地 生态 风险 21 年 间 总 体 呈 减 小 趋势 ,各 区 县 生态 风险 差异 显著 。 通 过 对 榆林 市 
土地 生态 风险 综合 评价 ， 可 为 优化 土地 资源 配置 、 保 护 生态 环境 提供 理论 建议 ， 实 现 区 域 经 济 生 态 和 谐 发 展 。 
关键 词 : 土地 ; 生态 风险 ; 生态 脆弱 ; 能 源 开发 ; 榆林 市 
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Abstract: Land use/cover change is a determinant factor in the analysis of ecosystem services for regional ecological security 
and maintenance. Ecological risks of land serves as a vital index in measuring sustainable utilization of land resources. The 
research area used in this study was Yulin City, a city in northern Shaanxi Province located in the Loess Plateau and Maowusu 
Desert border where the climate is arid and the land use structure is irrational. This region is one of the areas with the most 
severe water and soil erosion on the Loess Plateau, with ecologically-fragile environment, damage susceptibility and recovery 
difficulty. Thus based on the natural geographical conditions, development history of the city, basic land system characteristics 
and removed single factor evaluation uncertainty, this study established an evaluation index system of various sources risk in 
Yulin City. Percent precipitation anomaly index was used as indicator for the risk of drought, the risk of land use structure as 
the indicator for the effect of human activity on the ecological system, and wind/water erosion as the indicators of soil erosion. 
Then the four evaluation indicators were used to construct a comprehensive evaluation model for the qualitative and quantita- 
tive evaluation and analysis of the factors affecting ecological risk of land in Yulin City for the period 1990-2010. The main 
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conclusions of the study were as follows: (1) the distribution of precipitation was uneven, with percent precipitation anomaly 
index initially decreasing and later increasing in the period from 1990 to 2012 for Yulin City. The area with lower value of the 
index mainly covered the eastern part of the city. And the index increased from the south to the north. (2) There were large and 
significant changes in the differences in regional risks of land use structure. The land use structure risk first increased and then 
decreased in the course of the investigated 21-year period. (3) More severe soil water erosion occutred mainly in the eastern 
and the southern regions of the study area (driven mainly by topographic factors) with an initial mean increase and then a de- 
crease over the 1990-2012 period. (4) There was more severe wind erosion mainly in the northern region, with Yuyang District 
having the largest erosion and Suide and Mizhi Counties the smallest erosion, and with a decreasing trend in mean erosion for 
the investigated period. (5) The trend in ecological risk of land was decreased, and significant differences in ecological risk of 
land were observed among counties. The area of Yulin City with severe ecological risk in 1990 was mainly in the northwest, 
area with high ecological risk in 2001 was mainly in the central region and west of Yulin City (e.g., Dingbian County), area 
with moderate ecological risk was mainly in the southeast. The ecological security status in the study area was improved over- 
all, although the situation remained grim. This research put forth a theoretical proposal for optimizing the allocation of land 
resources and protecting the ecological environment in order to realize a harmonious development of regional economy and 
ecology. 
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近年 来 ， 随 着 人 口 增 长 和 社会 经 济 的 快速 发 展 ， 
人 类 的 社会 经 济 活动 导致 自然 资源 的 高 速 耗 减 和 环 


件 和 城镇 发 展 状况 ， 多 方面 综合 考虑 榆林 市 生态 环 
境 特征 ， 去 除 单一 要 素 评 价 的 不 确定 性 ， 建 立 多 种 


境 质量 下 降 , 由 此 产生 了 一 系列 生态 环境 问题 ， 严 
重 威胁 着 人 类 的 生存 和 发 展 。 人 类 该 如 何 维护 自身 
的 安全 成 为 人 们 关注 的 焦点 ， 生 态 风险 评价 也 越 来 
越 得 到 国内 外 的 关注 50“ 。 有 些 学 者 以 土地 利用 数 
据 为 基础 ,构建 生态 风险 指数 模型 ， 对 区 域 生态 风 
险 进行 了 综合 评价 "J 也 有 学 者 运用 景观 生态 学 
原理 对 区 域 生态 风险 进行 评价 "1 张 晓 媛 等 
基于 PSR 模型 构建 风险 源 、 风 险 受 体 和 风险 响应 评 
价 指标 体系 ,利用 综合 模糊 评价 法 对 三 峡 库 区 屏障 
带 重 庆 段 土地 利用 生态 风险 进行 综合 评价 ; 张 思 锋 
等 中 运用 生态 梯度 风险 评价 法 ,通过 定性 、 半 定 
量 、 定 量 3 个 层次 对 榆林 煤炭 开采 区 的 生态 风险 进 
行 评 价 。 

生态 风险 的 研究 最 早 主要 是 针对 特定 小 区 域 化 
学 污染 物 的 生态 风险 评价 0“ ， 随 着 经 济 发 展 、 人 
地 关系 矛盾 突出 , 土地 生态 风险 评价 逐渐 兴起 ,， 国 
内 学 者 结合 土地 利用 类 型 和 景观 因子 , 建立 了 数学 
评价 模型 和 基于 景观 尺度 的 生态 风险 框架 对 土地 生 
态 风 险 进行 评价 。 但 目前 对 于 土地 生态 风险 评价 指 
标 较为 单一 ， 多 局 限于 土地 利用 变化 ， 较 少 涉及 区 
域 自然 地 理 条 件 及 与 人 类 活动 相关 的 多 种 自然 、 社 
会 风险 ， 对 于 单一 生态 风险 评价 较为 成 熟 , 缺乏 在 
土地 生态 系统 结构 、 组 分 、 功 能 的 相互 作用 基础 上 
的 环境 效应 表征 ， 关 于 区 域 土地 生态 风险 综合 评价 
的 定量 化 和 空间 表达 仍 处 于 探索 阶段 。 

本 文 以 我 国 典 型 生态 脆弱 区 和 能 源 开发 区 一 一 
陕西 省 榆林 市 为 主要 研究 对 象 ， 依 据 其 自然 地 理 条 


风险 源 的 评价 指标 ,选取 表征 干旱 情况 的 降水 距 平 
百分率 指数 作为 干旱 风险 评价 指标 , 以 土地 利用 结 
构 风 险 作 为 土地 生态 安全 系统 自身 安全 性 的 评价 指 
标 , 在 榆林 市 北部 风沙 草滩 区 、 南 部 黄土 高 原 区 分 
别 选取 风力 侵蚀 和 水 力 侵蚀 作为 表征 土壤 侵蚀 的 生 
态 风险 指标 ,在 此 基础 上 构建 生态 风险 综合 评价 模 
型 ,定量 评价 1990 一 2012 年 榆林 市 土地 生态 风险 ， 
揭示 榆林 市 各 区 县 土地 生态 风险 的 空间 差异 及 其 时 


间 变 化 过 程 。 


1 研究 区 概况 

榆林 市 位 于 陕西 省 北端 ， 黄 河中 游 西 岸 ， 地 处 
107°14'51”"E~111°14’31"E, 36°48'58”"N~39°35'07”N, 
土地 总 面积 42 920.18 km”。 地势 西北 高 、 东 南 低 ， 地 
貌 类 型 以 古 长 城 为 界 ,北部 为 风沙 草滩 区 ， 南 部 为 
黄土 丘陵 沟 歼 区。 该 地 区 夏季 炎热 多 十 ,冬季 寒 准 
干燥 , 年 均 温 7~13 'C, 气温 日 较 差 较 大 , 年 均 降 水 
量 约 400 mm 左右 ， 季 节 性 差异 明显 。 榆林 市 是 我 国 
主要 的 能 源 基 地 ， 全 市 矿产 资源 极其 丰富 ， 被 誉 为 
“中 国 的 科威特 ”。2013 年 全 市 总 人 口 376.99 万 人 ， 
生产 总 值 (GDP)2 846.75 亿 元 , 工业 总 产值 1 943.56 
亿 元 , 占 GDP 的 68.2%。 


2 数据 来 源 与 研究 方法 
2.1 数据 来 源 

本 文 主要 数据 包括 : 1)1990 年 、2001 年 、2012 
年 土地 利用 数据 ,来 源 于 中 国 科 学 院 资源 环境 数据 
中 心 、 地 球 系 统 科学 数据 共享 平台 、 中 国 西部 环境 
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与 生态 科学 数据 中 心 以 及 野外 土地 调查 数据 ; 2) 气 
象 数据 主要 来 源 于 中 国 气象 科学 数据 共享 服务 网 ; 
3)DEM 数据 来 源 于 中 国 科 学 院 数据 应 用 环境 中 心 ; 
4) 道 路 交通 、 水 系 、 行 政 边 界 等 其 他 数据 ,来 自 于 
国家 基础 地 理 数 据 库 。 
2.2 ”研究 方法 
2.2.1 评价 指标 选取 

土地 生态 风险 评价 是 衡量 土地 资源 可 持续 利用 
的 重要 指标 , 产生 土地 生态 风险 的 因素 主要 有 自然 
因素 和 人 文 因素 , 榆林 市 自然 灾害 主要 有 旱灾 、 霜 
冻 、 暴 雨 、 大 风 、 冰 埠 等 ， 其 中 尤 以 旱灾 危害 最 为 
严重 ， 因 此 在 自然 因素 方面 , 选取 表征 干旱 情况 的 
降水 距 平 百分率 指数 作为 风险 评价 指标 。 人 文 因 素 
主要 有 人 口 城镇 化 、 能 源 开发 、 开 荒 造 田 、 森 林 砍 
伐 等 对 生态 系统 产生 的 间接 影响 ,土地 作为 人 类 活 
动 的 载体 ， 其 类 型 、 结 构 变 化 最 能 反映 人 类 活动 的 
类 型 及 强度 ,因此 选取 土地 结构 变化 这 一 指标 来 表 
征 人 类 社会 活动 对 生态 系统 产生 的 风险 。 此 外 , 由 
于 榆林 市 脆弱 的 自然 环境 有 太 人 类 活动 的 干扰 ， 使 其 
土壤 侵蚀 较为 严重 , 而 榆林 市 北部 为 风沙 草滩 区 ， 
主要 为 风力 侵 馈 ,南部 为 黄土 高 原 区 ,主要 为 水 力 
侵蚀 因此 选取 表征 土壤 侵蚀 的 风力 侵蚀 和 水 力 侵 
蚀 作为 生态 风险 评价 指标 。 
2.2.2 ”评价 指标 计算 方法 

1) 降 水 距 平 百分率 指数 

降水 距 平 百分率 指数 (P,) 反 映 了 某 一 时 段 的 
降水 量 与 同期 平均 降水 状态 的 偏离 程度 ， 其 计算 
模型 0” 为 : 


RP =x100% (1) 


式 中 :PP 为 某 一 时 段 的 降水 量 ; P 为 多 年 平均 降水 量 ， 
本 文中 P 取 1980 一 2010 年 30 年 的 平均 降水 量 。 降 
水 距 平 百分率 指数 是 表征 干旱 状况 的 指标 ,其 值 越 
大 表示 某 一 时 期 的 降水 量 高 于 平均 值 越 高 ,， 负 距 平 
表示 降水 量 低 于 平均 值 。 

基于 ArcGIS 采用 空间 插值 法 "得 到 连续 的 降 
水 量 空间 数据 ,为 了 保证 插值 的 准确 性 ,选取 榆林 
市 周边 宁夏 、 山 西 、 陕 西 、 甘 肃 、 内 蒙 的 58 个 站 点 
作为 插值 点 ， 生 成 1990 年 、2001 年 、2012 年 的 降 
水 数据 和 1980 一 2012 多 年 平均 降水 数据 , 将 3 年 的 
降水 图 层 分 别 与 1980 一 2012 多 年 平均 降水 图 层 进 
行 栅 格 计算 。 

2) 土 地 利用 结构 生态 风险 强度 指数 

土地 利用 结构 及 其 变化 在 一 定 程度 上 可 以 间接 
反映 土地 生态 风险 状态 ， 本 文 借鉴 已 有 的 研究 成 果 ， 
采用 的 土地 利用 结构 生态 风险 模型 外 如 下 : 

. B,xW 
BR (2) 

式 中 : ERI 为 土地 利用 结构 生态 风险 强度 指数 ; Bj 为 i 
种 地 类 的 面积 ; B 为 总 面积 ; WW 为 第 i 种 地 类 的 生态 
风险 权重 。 

借鉴 臧 淑 英 、 马 彩虹 等 扩 的 研究 成 果 采 用 层次 
分 析 法 对 各 地 类 进行 赋 权 ， 其 结果 见 表 1。 以 榆林 市 
土地 利用 数据 为 基础 ， 根 据 土地 利用 结构 风险 模型 
计算 各 区 县 土地 利用 结构 生态 风险 指数 。 


表 1 土地 利用 结构 风险 评价 中 不 同 土地 利用 类 型 的 权重 


Table 1 Weights of different land use types in the evaluation of ecological risk of land use structure 


土地 利用 类 型 耕地 林地 草地 水 域 建设 用 地 未 利用 地 

Land use type Cropland Forest land Grass land Water body Construction land Unused land 
态 风险 

生态 风险 权重 0.220 5 0.035 8 0.059 7 0.075 2 0.467 6 0.141 2 


Ecological risk weight 


3) 土 壤 水 力 侵蚀 量 
土壤 水 力 侵蚀 估算 模型 中 ,目前 应 用 最 广泛 的 
两 个 模型 为 土壤 流失 通用 方程 (USLE) 和 RUSLE 模 
型 ， 后 者 是 在 前 者 修正 后 的 基础 上 得 到 的 ， 其 对 
USLE 中 的 一 些 因子 做 了 进一步 的 分 解 ， 估 算 更 为 
科学 、 精 确 ， 因 此 本 文采 用 RUSLE 模型 对 土壤 水 力 
侵蚀 量 进 行 估算 。 其 模型 基本 公式 为 ?1 
A, =RxKxCxLSxP (3) 
式 中 : 4 为 研究 区 域 土壤 实际 的 侵蚀 量 , 单位 为 
thm “a '; R 为 降水 -径流 因子 , 用 多 年 平均 年 降水 


侵蚀 力 指 数 表 示 ,， 单位 为 MJjmm-hm ha!; KK 为 
土壤 可 侵蚀 性 因子 ,单位 为 th.MJ mm ; C 为 覆盖 
与 管理 因子 ; LS 为 坡 长 坡度 因子 ; P 为 土壤 保持 措施 
因子 。C、LS、PP 均 为 无 量 纲 因 子 。 

4) 土 壤 风 力 侵蚀 量 

针对 不 同时 空 尺度 的 土壤 风力 侵 刨 情况， 出现 
了 WEQ 模型 WEPS 模型 \ 日 风沙 侵蚀 量 模型 等 , 但 
大 多 数 模型 都 难以 估算 长 时 间 的 累积 侵蚀 量 。 董 治 
宝 等 中 根据 我 国 的 实际 情况 , 建立 了 年 尺度 的 区 
域 土壤 风力 侵蚀 估算 模型 ,其 表达 式 为 : 
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o= 咱 加 1.0413+0.04410+0.002102_0.000103 


bXy 


人 
(82405) CRspr? /Hq? | dedydt 
式 中 : 2 为 土壤 风力 侵蚀 量 , t; VV 为 风速 , ms '; 及 为 
空气 相对 湿度 , %; VCR 为 植被 覆盖 度 , %; SDR 为 人 
为 地 表 结 构 破 损 率 , %; a 为 颗粒 平均 粒 径 , mm; FF 为 
土 体 硬度 , N.cm-; 0 为 坡度 , °; x， y 表示 距离 参照 点 
的 距离 , km; t 为 时 间 , s。 

5) 土 地 生态 风险 综合 指数 

在 参考 大 量 文献 以 及 结合 榆林 市 实际 情况 的 基 


础 上 ”4， 通 过 构建 生态 风险 综合 指数 对 榆林 市 生 
态 风 险 进行 评价 ， 其 模型 如 下 : 


E= > (Wx8) (5) 
i=] 


式 中 : 为 土地 生态 风险 综合 指数 ; 友 ; 为 单项 风险 权重 ; 
Pi 为 单项 生态 风险 评价 值 。 其 中 到 ;是 各 评价 指标 的 权 
重 值 ， 本 文采 用 定性 和 定量 相 结 合 的 层次 分 析 法 , 结合 
Yaahp 软件 确定 各 生态 风险 指标 的 权重 ( 表 2)。 根 据 土 
地 生态 风险 综合 指数 模型 计算 得 到 区 域 的 土地 生态 风 
伶 综合 指数 ,结合 已 有 研究 ,将 土地 生态 风险 划分 为 5 
个 等 级 ,其 划分 标准 如 表 3 所 示 。 


表 2 土地 生态 风险 综合 评价 中 各 分 类 生态 风险 的 权重 


Table 2 Weights of different kinds of ecological risks in the comprehensive evaluation of land ecological risk 


风险 干旱 风险 土地 利用 结构 风险 水 土 流失 风险 风沙 侵蚀 风险 
Ecological risk Drought risk Land use structure risk Water erosion risk Wind erosion risk 
T™ 权重 Weight 0.122 2 0.423 0 0.227 4 0.227 4 
表 3 土地 生态 风险 等 级 划分 标准 
Table 3 Classification standards of land ecological risk 
\ 等 级 Class | ll 咱 IV V 
名 生态 风险 指数 Ecological risk index <0.25 0.25~0.35 0.35~0.45 0.45~0.55 >0.55 
风险 状态 Risk state 很 小 Tiny 较 小 Less 稍 小 Small 较 大 Larger 很 大 Great 


3 结果 与 分 析 

干旱 风险 评价 

通过 计算 得 到 1990 年 、2001 年 和 2012 年 的 降 
水 距 平 百分率 指数 分 布 图 (图 1), 为 了 进一步 定量 化 
分 析 ， 以 县 域 为 行政 单元 ， 进 行 分 区 统计 ， 得 到 各 
区 县 的 降水 距 平 百 分 率 指数 ( 表 4)。 由 表 4 和 图 1 可 
知 : 榆林 市 1990 年 的 降水 距 平 百分率 指数 低 值 区 主 
要 分 布 在 东北 部 ， 高 值 区 主要 分 布 在 南部 ， 其 中 低 
值 中 心 达 -0.15%， 高 值 中 心 为 29.52%， 大 致 呈现 由 


3.1 


1990 2001 
| | High: 29.51 国 High: 48.31 
Low: —0.15 Low: —16.96 


北向 南 的 递增 趋势 。 各 区 县 中 ， 府 谷 县 、 神 木 县 距 
平 值 较 低 ， 而 清润 县 最 高 。 整 个 年 份 大 部 分 地 区 均 
为 正 值 ， 表 明 其 降水 量 在 均值 以 上 , 但 高 值 和 低 值 
差异 较 大 ， 说 明 榆 林 市 各 区 县 降水 分 布 不 均衡 ， 表 
现 为 南 多 北 少 。2001 年 的 降水 距 平 百分率 指数 低 值 
区 主要 分 布 在 东部 ,高 值 区 主要 分 布 在 中 部 和 西部 ， 
大 致 呈现 自 东 向 西 递增 的 趋势 ， 其 中 距 平 指数 最 低 
的 是 府 谷 县 ， 最 高 值 为 榆 阳 区 。 该 年 份 低 值 中 心 达 
-16.96%， 相 上 比 1990 年 明显 降低 , 说 明 2001 年 降水 
量 减少 , 干旱 程度 加 剧 。 但 高 值 中 心 为 48.31%， 高 
N 


人 


2012 


High: 68.35 
Low: —9.83 


100 km 


1 
Distribution of precipitation anomaly percentage index (%) of Yulin City in 1990, 2001 and 2012 


1990 年 、2001 年 和 2012 年 榆林 市 降水 距 平 百 分 率 指数 (%) 分 布 图 
Fig. 1 
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值 与 低 值 差距 增 大 , 说明 降 水 分 布 不 均衡 加 剧 ， 可 
能 同时 发 生 旱 涝 ; 2012 年 的 降水 距 平 百分率 指数 低 
值 区 主要 分 布 在 东南 部 ,高 值 区 主要 分 布 在 东北 部 ， 
大 致 呈现 由 南 向 北 递增 的 趋势 ， 其 中 距 平 指数 较 低 
的 为 清润 县 ， 较 高 的 为 府 谷 县 。 与 1990 年 、2001 年 
相 比 ， 其 降水 距 平 指数 值 整体 增加 ,说 明 2012 年 的 
降水 量 明 显 增加 ,高 值 与 低 值 差异 明显 ， 降 水 分 布 
仍然 不 均衡 。 

1990 年 、2001 年 和 2012 年 榆林 市 各 区 县 

降水 距 平 百分率 指数 


Table 4 Precipitation anomaly percentage index of each 
county (district) of Yulin City in 1990, 2001 and 2012 % 


表 4 


区 县 County (district) 1990 2001 2012 
米 脂 县 Mizhi County 16.23 8.35 23.87 
绥 德 县 Suide County 24.57 7.95 13.57 
靖 边 县 Jingbian County 21.73 6.71 25.36 
榆 阳 区 Yuyang District 3.52 35.13 44.98 
佳 县 Jiaxian County 11.09 7.03 27.57 
府 谷 县 Fugu County 2 —8.24 55.97 
神木 县 Shenmu County .62 2.44 44.84 
定 边 县 Dingbian County 23.72 26.49 24.38 
子 州 县 Zizhou County 21.63 4.43 20.53 
清润 县 Qingjian County 28.10 2.49 11.27 
吴 堡 县 Wubao County 24.13 2.47 13.05 
横山 县 Hengshan County 13.58 23.81 35.03 
均值 Mean 15.93 14.09 28.37 
表 5 


3.2 土地 利用 结构 风险 评价 

由 表 5 可 知 : 1990 年 土地 利用 结构 风险 指数 最 
大 的 是 府 谷 县 ， 最 小 的 是 绥 德 县 。 土 地 利用 生态 风险 
较 大 的 县 域 主要 分 布 在 北部 ， 其 中 包括 府 谷 县 、 神 木 
县 、 榆 阳 区 ， 而 位 于 东南 部 的 佳 县 、 米 脂 县 、 吴 堡 县 、 
绥 德 县 和 子 洲 县 土地 利用 结构 生态 风险 较 小 。 产 生 这 
一 分 布 的 原因 主要 与 北部 各 县 耕地 、 未 利用 地 所 占 比 
重 较 大 , 而 生态 性 用 地 林地 、 草 地 、 水 域 所 占 比 重 较 
小 有 关 。2001 年 土地 利用 结构 生态 风险 整体 分 布 发 
生 明 显 变化 , 土地 利用 结构 风险 指数 最 大 的 是 米 脂 
县 , 最 小 的 是 神木 县 ,， 生态 风险 较 1990 年 明显 加 剧 ， 
这 主要 由 于 城市 迅速 扩张 ,建设 用 地 上 比重 增加 ， 其 中 
子 洲 县 、 米 脂 县 、 绥 德 县 、 佳 县 和 吴 堡 县 生态 风险 显 
著 加 剧 ， 主 要 与 其 各 地 类 转移 , 且 耕 地 面积 所 占 比重 
明显 增加 有 关 。2012 年 土地 利用 结构 生态 风险 较 大 
的 县 域 主要 分 布 在 南部 ， 其 中 指数 最 大 的 是 吴 堡 县 ， 
最 小 的 是 榆 阳 区 。 从 土地 利用 结构 风险 指数 变化 来 看 ， 
1990 一 2001 年 , 土地 利用 结构 风险 指数 减少 的 有 府 
谷 县 、 靖 边 县 、 神 木 县 、 榆 阳 区 ， 其 余 各 县 风险 指 
数 均 有 所 增加 ; 2001 一 2012 年 , 土地 利用 结构 风险 
旨 数 增加 的 有 府 谷 县 、 米 脂 县 、 绥 德 县 、 清 润 县 、 
吴 堡 县 ， 其 余 各 县 的 风险 指数 均 有 所 减少 ; 从 总 体 
来 看 , 1990 一 2012 年 ， 榆 阳 区 、 神 木 县 、 靖 边 县 的 土 
地 利用 结构 风险 指数 一 直 处 于 减少 趋势 , 米 脂 县 、 
清润 县 、 绥 德 县 、 吴 堡 县 的 土地 利用 结构 风险 指数 
一 直 处 于 增加 趋势 。 


1990 年 、2001 年 和 2012 年 榆林 市 土地 利用 结构 风险 指数 及 其 变化 


Table S$ Risk indexes and changes of land use structure in different counties (district) in Yulin City in 1990, 2001 and 2012 


土地 利用 结构 风险 指数 土地 利用 结构 风险 指数 变化 
区 县 County (district) Land use Structure risk index Change of land use Structure risk index 
1990 2001 2012 1990 一 2001 2001 一 2012 1990 一 2012 

定 边 县 Dingbian County 0.122 6 0.139 5 0.1317 0.016 9 -0.0079 0.009 0 
府 谷 县 Fugu County 0.158 7 0.117 1 0.121 3 —0.041 6 0.004 2 —0.037 4 
横山 县 Hengshan County 0.1161 0.1463 0.116 1 0.030 2 -0.030 2 0.000 1 
佳 县 Jiaxian County 0.102 4 0.154 5 0.150 7 0.052 1 -0.003 8 0.048 3 
靖 边 县 Jingbian County 0.132 0 0.127 1 0.110 4 —0.004 8 -0.016 7 -0.021 6 
米 脂 县 Mizhi County 0.097 2 0.159 0 0.160 8 0.061 8 0.0019 0.063 7 
清润 县 Qingjian County 0.117 3 0.1266 0.141 8 0.009 4 0.015 2 0.024 5 
神木 县 Shenmu County 0.146 5 0.114 6 0.092 6 -0.0319 -0.022 0 -0.053 9 
绥 德 县 Suide County 0.089 0 0.149 2 0.162 8 0.060 2 0.013 6 0.073 8 
吴 堡 县 Wubao County 0.105 5 0.148 1 0.173 3 0.042 7 0.025 2 0.067 8 
榆 阳 区 Yuyang District 0.141 6 0.123 9 0.086 5 -0.017 7 —0.0374 —0.055 1 
子 洲 县 Zizhou County 0.105 9 0.145 7 0.133 2 0.039 8 -0.012 5 0.027 3 


3.3 ”土壤 水 力 侵蚀 评价 
根据 土壤 水 力 侵蚀 模型 将 所 获取 的 各 因子 图 层 
进行 栅 格 二 加 计算 ， 得 到 榆林 市 1990 年 、2001 年 和 


2012 年 的 土壤 水 力 侵蚀 强度 图 (图 2)。 从 图 2 可 知 ， 
侵蚀 严重 的 地 区 主要 分 布 在 榆林 市 东北 部 及 中 南部 ， 
东北 部 处 于 水 蚀 - 风 蚀 交 错 区 ， 中 南部 地 形 复杂 ， 坡 
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度 大 , 土壤 侵蚀 严重 。 从 时 间 变 化 看 , 1990 年 榆林 市 
东南 部 地 区 土壤 水 力 侵蚀 较 严 重 ， 其 最 高 侵蚀 量 达 
1 692.49 t:hm a1, 平均 侵蚀 量 约 147.06 thm a， 
其 中 侵蚀 量 最 大 的 为 清润 县 , 侵蚀 量 最 小 的 为 定 边 
县 ; 2001 年 榆林 市 东部 及 南部 土壤 水 力 侵蚀 均 较 严 
重 ， 其 最 高 侵蚀 量 达 3 711.99 thm-”.a 1， 平均 侵蚀 
量 约 为 296.65 thm a '， 其 中 府 谷 县 侵蚀 强度 最 大 ， 
榆 阳 区 侵蚀 量 最 小 ; 2012 年 榆林 市 南部 土壤 水 力 侵 
蚀 较 严重 ， 其 最 高 侵蚀 量 达 3 076.8 thm a!, 平均 
侵蚀 量 约 为 181.10 thm*a-!， 其 中 清润 县 侵蚀 强度 


1990 2001 


High: 1 692.49 
Low: 0 


‘High: 3 711.99 


四 Low: 0 


100 km 


最 大 ， 榆 阳 区 侵蚀 强度 最 小 。 将 1990 年 、2001 年 和 
2012 年 土壤 水 力 侵蚀 强度 进行 比较 ,侵蚀 强度 2001 
年 >2012 年 >1990 年 。1990 一 2001 年 土壤 侵蚀 强度 
明显 加 剧 ， 其 中 府 谷 县 、 靖 边 县 、 定 边 县 最 为 明显 ; 
到 2012 年 , 各 县 土壤 水 力 侵蚀 强度 有 所 减缓 ， 变 化 
最 为 显著 的 是 府 谷 县 、 神 木 县 、 定 边 县 ; 总 体 来 看 ， 
1990 一 2012 年 神木 县 、 府 谷 县 、 吴 堡 县 、 佳 县 的 土 
塘 水 力 侵蚀 强度 有 所 缓和 ， 其 余 各 县 的 土壤 水 力 侵 
蚀 强 度 均 呈现 不 同 程度 的 加 剧 ， 其 中 清润 县 、 靖 边 
县 和 子 洲 县 最 为 显著 。 


N 
2012 人 


High: 3 076.80 


图 Low 0 


2 1990 年 、2001 年 和 2012 年 榆林 市 土壤 水 力 侵蚀 强度 (thm-?.a-0 分 布 图 
Fig. 2 Distribution of soil erosion intensities (thm ”a !) of Yulin City in 1990, 2001 and 2012 


3.4 土壤 风力 侵蚀 评价 

根据 土壤 风力 侵蚀 模型 将 所 获取 的 各 因子 图 层 
进行 栅 格 盖 加 计算 , 最 终 得 到 榆林 市 1990 年 、2001 
年 和 2012 年 的 土壤 风力 侵蚀 强度 图 (图 3)。 如 图 3 
所 示 ，1990 年 榆林 市 西部 及 北部 地 区 风力 侵蚀 较 严 
重 ， 其 最 高 侵蚀 量 达 128.83 tkm“, 平均 侵蚀 量 约 为 
22.91 tkm“， 其 中 榆 阳 区 和 定 边 县 侵蚀 量 最 大 ， 绥 德 
县 和 清润 县 侵蚀 量 最 小 ; 2001 年 榆林 市 北部 地 区 风力 


1990 2001 


High: 128.833 1 
Low: 0.004 7 


‘High: 130.398 0 
国 1 0.009 8 


0 


100 km 


侵蚀 仍 较为 严重 ， 其 最 高 侵蚀 量 达 130.398 tkm , 平 
均 侵 蚀 量 约 为 2.82 t-km“, 榆 阳 区 和 神木 县 侵蚀 量 最 
大 , 绥 德 县 和 米 脂 县 侵蚀 量 最 小 ; 2012 年 榆林 市 北部 
地 区 风力 侵蚀 较 严 重 ， 其 最 高 侵蚀 量 达 52.12 tkm 
平均 侵蚀 量 约 为 0.70 tkm”， 其 中 榆 阳 区 侵蚀 量 最 大 ， 
其 次 为 神木 县 ,， 米 脂 县 和 绥 德 县 侵蚀 量 最 小 。 将 1990 
年 、2001 年 和 2012 年 土壤 风力 侵蚀 强度 进行 比较 , 侵 
蚀 强 度 1990 年 >2001 年 >2012 年 。 从 侵蚀 强度 分 布 图 


N 
2012 人 


‘High: 52.118 2 


ml Low: 0.003 6 


3 1990 年 、2001 年 和 2012 年 榆林 市 土壤 风力 侵蚀 强度 (tkm- ”分布 图 
Fig. 3 Distribution of sand erosion intensities (tkm- *) of Yulin City in1990, 2001and 2012 


http:/www.ecoagri.ac.cn 


662 


中 国生 态 农业 学 报 2017 Ch 


inaXiv 合 人 仍 汪 前 


来 看 ,侵蚀 严重 的 地 区 主要 分 布 在 榆林 市 北部 地 区 ， 
主要 由 于 北部 位 于 风沙 草滩 区 ,植被 覆盖 率 低 。 从 
侵蚀 量变 化 来 看 ，1990 一 2012 年 整体 呈 减 少 趋势 ， 
减少 量 最 多 为 榆 阳 区 ， 其 次 为 定 边 县 、 横 山 县 ， 其 中 
1990 一 2001 年 相 较 于 2001 一 2012 年 侵蚀 强度 变化 
明显 ， 其 年 均 变 化 率 约 为 2001 一 2012 年 的 1.2 倍 。 
3.5 土地 生态 风险 综合 评价 

由 表 6 可 知 : 1990 年 土地 生态 风险 综合 指数 最 
大 的 是 府 谷 县 ， 最 小 的 是 绥 德 县 ， 两 者 差异 明显 ， 
最 大 值 约 为 最 小 值 的 2.8 倍 ; 2001 年 土地 生态 风险 综 
合 指数 最 大 的 是 佳 县 ， 最 小 的 是 神木 县 ， 最 大 值 约 
为 最 小 值 的 2 倍 ， 有 所 减 小 ; 2012 年 土地 生态 风险 综 
合 指数 最 大 的 是 清润 县 ， 最 小 的 是 神木 县 ， 其 最 大 
值 约 为 最 小 值 的 1.9 倍 , 较 1990 年 和 2001 年 差异 有 
所 减 小 。 从 各 区 县 变化 来 看 ，1990 一 2001 年 ， 靖 边 
县 、 榆 阳 区 、 府 谷 县 、 神 木 县 和 清润 县 的 土地 生态 
风险 综合 指数 减少 ,说 明 其 生态 风险 有 所 改善 ， 其 


中 减少 最 多 的 是 府 谷 县 ， 其 年 递减 率 约 为 4.48%6; 
其 余 各 县 的 土地 生态 风险 综合 指数 均 有 所 增加 ， 说 
明 其 生态 风险 加 剧 ， 其 中 增加 最 多 的 是 米 脂 县 ， 其 
年 递增 率 约 为 11.99%。2001 一 2012 年 , 府 谷 县 、 清 
润 县 和 吴 堡 县 的 土地 生态 风险 综合 指数 增加 ， 其 中 
增加 最 多 的 是 清润 县 ， 其 年 递增 率 约 为 3.73%， 其 
余 各 县 的 土地 生态 风险 综合 指数 均 有 所 减少 ， 其 中 
减少 最 多 的 是 横山 县 ， 其 年 递减 率 约 为 3.74%。 从 总 
体 来 看 , 1990 一 2012 年 ， 米 脂 县 、 绥 德 县 、 佳 县 、 子 
洲 县 、 吴 堡 县 的 土地 生态 风险 综合 指数 增加 ， 其 中 
增加 最 多 的 是 绥 德 县 ， 其 年 递增 率 约 为 4.78%， 其 
余 各 县 的 土地 生态 风险 综合 指数 减少 ， 其 中 减少 最 
多 的 是 神木 县 ， 其 年 递减 率 约 为 2.50%。 

由 图 4 可 知 ，1990 一 2012 年 榆林 市 土地 生态 风 
险 等 级 发 生 了 较 大 变化 。1990 年 榆林 市 各 区 县 中 前 
3 等 级 的 县 域 共 有 5 个 , 约 42%， 主 要 位 于 东南 部 地 
区 , 土地 生态 风险 较 小 ， 生 态 环境 较为 安全 。 而 位 于 


表 6 1990 年 、2001 年 和 2012 年 榆林 市 土地 生态 风险 综合 指数 及 其 变化 
Table6 Land comprehensive risk indexes and changes in different counties (district) in Yulin City in 1990, 2001 and 2012 


土地 生态 风险 综合 指数 
区 县 County (district) 


Land comprehensive risk index 


土地 生态 风险 综合 指数 变化 


Change of land comprehensive risk index 


1990 2001 2012 1990 一 2001 2001 一 2012 1990 一 2012 

米 脂 县 Mizhi County 0.262 8 0.609 2 0.478 9 0.346 4 -0.130 3 0.216 2 
绥 德 县 Suide County 0.239 1 0.495 7 0.490 5 0.256 7 -0.005 2 0.251 4 
靖 边 县 Jingbian County 0.516 7 0.398 3 0.290 8 -0.118 5 -0.107 5 -0.226 0 
榆 阳 区 Yuyang District 0.561 0 0.438 2 0.3196 -0.122 8 -0.118 6 —0.241 4 
佳 县 Jiaxian 0.316 8 0.639 0 0.447 3 0.322 1 -0.1916 0.1305 
府 谷 县 Fugu County 0.659 9 0.334 9 0.444 3 -0.325 1 0.109 5 -0.215 6 
神木 县 Shenmu County 0.600 2 0.323 3 0.269 7 -0.276 9 -0.053 6 -0.330 5 
定 边 县 Dingbian County 0.492 7 0.496 5 0.340 1 0.003 8 -0.156 3 -0.152 6 
子 洲 县 Zizhou County 0.406 4 0.530 5 0.427 7 0.124 1 -0.102 8 0.021 3 
清润 县 Qingjian County 0.535 3 0.359 5 0.506 9 -0.175 8 0.147 4 —0.028 4 
吴 堡 县 Wubao County 0.400 6 0.477 5 0.498 5 0.0769 0.021 1 0.0979 
横山 县 Hengshan County 0.466 2 0.593 5 0.349 7 0.127 3 -0.243 9 -0.116 5 
N 
1990 2001 2012 人 

I 级 Class I 

[7 [级 Class 了 

亚 级 Class 亚 

区 级 ClassIV 

国 量 V 级 Class V 


图 4 


1990 年 、2001 年 和 2012 年 榆林 市 生态 风险 等 级 分 布 
Fig.4 Distribution of ecological risk grade of Yulin City in 1990, 2001 and 2012 
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西部 毛乌素 沙 地 区 生态 风险 较 大 ， 该 区 域 水 土 流 
严重 ， 植 被 覆盖 度 较 低 ， 生 态 环境 遭受 严重 破坏 ; 
从 2001 年 生态 风险 等 级 分 布 图 看 , 各 区 县 中 没有 属 
于 | 级 生态 风险 等 级 的 县 ，IV 级 、V 级 的 县 域 共 有 7 
个 , 约 占 58%， 主 要 位 于 榆林 市 东南 部 以 及 西北 的 
定 边 县 ， 生 态 环境 仍然 处 于 不 安全 状态 ; 2012 年 , 全 
市 各 区 县 中 没有 属于 | 级 、V 级 生态 风险 等 级 的 县 ， 
其 中 前 三 等 级 的 县 域 个 数 共 有 8 个 , 约 占 67%， 较 
1990 年 、2001 年 相 比 ， 生 态 风险 有 所 减 小 ,但 是 仍 
然 有 4 个 县 域 生态 风险 较 大 ， 生 态 风险 依然 严峻 。 
4 结论 

本 文选 取 干 旱 、 土 地 利用 结构 、 土 壤 水 力 侵蚀 、 
风力 侵蚀 4 个 风险 指标 ， 构 建 区 域 生 态 风险 综合 评 
价 模型 ， 对 榆林 市 1990 一 2012 年 土地 生态 风险 进行 
探索 性 评价 ,研究 结果 如 下 : 

1) 榆 林 市 各 生态 风险 要 素 时 空 差异 较 大 。1990 一 
2012 年 间 榆 林 市 的 降水 距 平 百分率 指数 均值 先 减 少 后 
增加 ， 大致 呈 现 由 南 向 北 递增 的 趋势 ， 其 低 值 区 主要 
分 布 在 偏 东 部 地 区 , 高 值 区 与 低 值 区 差异 明显 , 降水 
分 布 不 均衡 ; 土地 利用 结构 生态 风险 区 域 差异 较 大 和 且 
变化 明显 ， 呈 先 增 后 减 的 趋势 ,其 中 米 脂 县 、 清 润 县 、 
绥 德 县 、 吴 堡 县 的 土地 利用 结构 风险 指数 一 直 处 于 增 
加 趋势 ， 需 要 引起 重视 ， 而 榆 阳 区 一 直 处 于 减少 趋 
势 ， 情况 有 所 改善 ; 土壤 水 力 侵蚀 较 剧烈 的 区 域 主 
要 分 布 在 东部 及 南部 地 区 ,这 主要 与 该 区 域 地 形 因 
素 有 关 ， 其 均值 呈 先 增 后 减 的 趋势 ， 其 中 神木 县 、 府 
谷 县 、 佳 县 和 吴 堡 县 的 土壤 水 力 侵蚀 强度 有 所 缓和 ， 
而 清润 县 、 靖 边 县 和 子 洲 县 呈现 显著 加 剧 现象 ; 风 
力 侵蚀 较 严 重 的 区 域 主要 分 布 在 北部 地 区 ， 侵 蚀 量 
最 大 的 是 榆 阳 区 ， 最 小 的 为 绥 德 县 、 米 脂 县 ， 其 均值 
一 直 呈 递减 趋势 ， 榆 阳 区 减少 量 最 多 ， 其 次 为 定 边 
县 、 横 山 县 。 

2) 榆 林 市 土地 生态 风险 等 级 整体 变化 较 大 且 区 
域 差异 明显 。1990 年 生态 风险 等 级 属于 IV、V 级 的 
县 域 有 7 个 ,生态 风险 较 大 ， 其 中 风险 最 大 的 是 府 
谷 县 ; 2001 年 生态 风险 等 级 属于 IV、V 级 的 县 域 仍 
然 有 7 个 ,生态 风险 较 大 ， 生 态 环境 处 于 安全 状态 ， 
其 中 风险 最 大 的 是 佳 县 ; 2012 年 生态 风险 有 所 减缓 ， 
但 是 生态 风险 等 级 属于 IV、V 级 的 县 域 仍 有 4 个 ， 生 
态 安 全 依然 严峻 。 从 区 域 差异 来 看 ，1990 年 榆林 市 
生态 风险 较 大 的 区 域 主要 位 于 西北 部 地 区 ,2001 年 
生态 风险 较 大 的 区 域 主要 位 于 榆林 市 中 部 地 区 和 西 
边 的 定 边 县 ，2012 年 生态 风险 较 大 的 地 区 主要 分 布 


在 东南 部 ， 生 态 安 全 状态 整体 有 所 改善 ,但 形势 仍 
旧 严 峻 ， 需 要 引起 注意 。 

榆林 市 属于 干旱 半 干 旱 区 、 典 型 的 农 牧 交 错 带 
和 能 源 开 发 基地 ， 生 态 环境 脆弱 ， 人 地 矛盾 突出 
受 能 源 经 济 产业 化 和 生态 恢复 过 程 的 双重 影响 ,土地 
生态 安全 问题 越 来 越 突出 。 本 文 从 自然 灾害 、 人 为 活 
动 和 生态 背景 出 发 , 构建 了 区 域 尺度 的 综合 生态 风险 
评价 指标 体系 和 估算 模型 ， 定量 评价 了 1990 一 2012 年 
榆林 市 土地 生态 风险 的 时 空 变化 差异 ， 有 针对 性 地 分 
析 了 其 动态 变化 趋势 。 结 果 表明 随 着 榆林 市 能 源 的 大 
量 开 采 和 煤矿 产业 的 莲 勃 发 展 ， 榆林 市 生态 环境 遭受 
了 严重 的 破坏 ， 提 高 生态 环境 意识 , 加 强生 态 文 明 建 
设 至 关 重要 。 下 一 步 研究 将 围绕 土地 生态 风险 的 发 
生机 制 建立 合理 的 土地 生态 风险 评价 指标 体系 , 细 
化 区 域 研究 尺度 ， 提 高 空间 表达 的 精度 ， 运 用 主观 
和 客观 赋 权 相 结 合 合理 确定 指标 权重 ， 从 评价 指标 
体系 的 构建 、 区 域 研究 尺度 的 细 化 等 方面 提高 土地 
生态 风险 综合 评价 的 科学 性 , 为 合理 进行 土地 资源 
管理 规划 提供 决策 。 
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